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FOREWORD 


This  is  the  ninth  report  in  the  Data  Reduction  Program  Documentation  series. 

It  is  dated  according  to  the  date  of  completion  of  the  documentation.  No  implication 
is  made  that  this  program  will  not  subsequently  be  modified,  amended,  or  superseded, 
on  the  contrary,  the  history  radar  data  processing  is  one  of  continuous  evolution  o 
techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has  been  reached.  The 
PA-229  series  is  being  published  for  the  convenience  of  interested  parties,  and  Lincoln 
assumes  no  responsibility  for  the  correctness  of  the  information  presented,  nor  for 
its  currency. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles R. 
Berndtson  of  Lincoln,  and  of  D.  Nessman  and  R.  French  of  Philco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Berndtson. 

The  principal  contributor  to  this  report  was  G.  L.  Shapiro  (Philco-Ford).  Due 
to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being,  the 
editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  --  mathematicians 
or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and  writing  of 

the  programs. 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 


(The  units  given  for  certain  quantities  are  the  units  commonly  usert  for  those 
quantities,  unless  otherwise  noted. ) 


ADT 

ALCOR 

ALTAIR 

Alt 

APS 

ARS 

Avg 

Az 


ALCOR  Data  Tape 

ARPA-Lincoln  C-band  Observables  Radar 

ARPA  Long-Range  Tracking  and  Instrumentation  Radar 

Altitude  (km) 

Average  Fulse  Shape 
ALTAIR  Recording  System 
Average,  Averaging 
Azimuth  (deg) 


c 

CADJ 
C -band 

DBLT 

El 

EOF 


Speed  of  Light 

Adjusted  Calibration  Constant  (db) 

ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 


Wide  Band  Pulse  Doublet 


Elevation  (deg) 
End  of  File 


GMT 


Greenwich  Mean  Time 


h 

Hz 

IF 

in 

LC 

lsb 


Hours 

Hertz 


Intermediate  Frequency 
Inches 

Left  Circular  Polarization 
Least  Significant  Bit 


min 

NB 

nrtpod 

POD 

Phase 

PRF 

PRI 

pps 

pts 


Minutes 
Narrow  Band 

Non -real  Time  Precision  Orbit  Determination  Program 

Project  PRESS  Operation  and  Data  Summary  Report 
Presented  in  deg 

Pulse  Repetition  Frequency  (pps) 

Pulse  Repetition  Interval  (s) 

Pulses  per  second 
Points 


v 


rad 

RC 

RCS 

RF 


Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 
Radio  Frequency 


SDBLT 

S/N 

T 

TAL 

UHF 

V 


V 

w 

VHF 

WB 

WBS 

WTR 

e 

x 


Seconds 

Standard  Deviation  of  Wake  Velocity 

Wide  Band  Slaved  Pulse  Doublet 
Signal -to -noise  Ratio 

Time 

Time  After  Launch  (s) 

ALTAIR  Frequency;  415  MHz 
Velocity 

Doppler  Velocity 

Mean  Wake  Velocity 

ALTAIR  Frequency;  155.5  MHz 

Wide  Band 

Wide  Band  Slaved 

Western  Test  Range 

Total  Off-axis  Angle  (deg) 
Wavelength 

Denotes  Multiplication 


FLOW  DIAGRAM  SYMBOLS 


PROCESS,  ANNOTATION 


TERMINATOR 


SUBROUTINE:  where  NAME  Is  the  entry 
call  Into  the  subroutine 


CONNECTOR:  where  P  speclfles  a  page  In  the 

flow  diagram,  and  E  designates 
a  statement  number  In  the  program 
listing  or  a  reference  point  in  the 
flow  diagram 

CONNECTOR:  where  X  implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


PUNCHED  CARD 


DISK 
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ALC10 


I. 


PURPOSE  and  UTILIZATION 


A. 


B. 


Source  of  Data 


ALCOR 


Data  Input 

ALCOR  Data  Tape  (ADT) 
description 


(jf/k-Ji-'OSK.  dA-fc-**  ^aJirurt^rOrOMiit  om£,. 

ALC10  is  designed  to  ootainJRCS-data  on  hard  body  targets  (NB-apd  WS) 
and  on  multiple  radar  scattering  centers  {WS-).  It  produces  plots  vs^FAL  and -Alt  of 


the  peak  RCS  from  a  set  of  gates  or  of  the  RCS  in  a  selected  gate 
D.  Output 


K 


\ 


1.  A  listing  of  pertinent  identification  da*-". 

2.  Plots  vs  TAL  and  Alt  of  peak  RCS  or  of  RCS  in  a  selected 
range  gate. 
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II. 


DESCRIPTION 


ALC10  computes  average  RCS  as  a  function  of  range  gate  and  time.  The 
2 

program  averages  in  m  and  then  converts  to  dbsm  for  printouts  and  plots. 

The  averaging  interval  (A VINT)  is  input  in  seconds,  and  must  be  either 
0.01,0.02,0.04,0.05,0.08,  or  0.1.  Based  on  AVINT  and  the  primary  pulse  PRF, 
the  number  of  pulses  averaged  is  shown  in  Table  I. 

Only  single  pulse  RCS  data  are  computed  if  the  range  offset  (XINOFF) 
is  not  zero. 


RCS  is  computed  as  follows: 


LC  RCS  =  XATBL(N)  +  40  log  R  +  XPPAGC  +  CONLC  -  POWERT 
RC  RCS  =  XATBL(M)  +  40  log  R  +  XOPAGC  +  CONRC  -  POWERT 

where 

XATBL(N)  is  obtained  by  indexing  the  LC  calibration  table  with  th^ 
LC  amplitude  values  obtained  in  the  ADT  data  record. 

XATBL(M)  is  obtained  by  indexing  the  RC  calibration  table  with  the 
RC  amplitude  values  obtained  in  the  ADT  data  record. 1 

XPPAGC  is  otal  LC  attenuation  (db)2 

XOPAGC  is  total  RC  attenuation  (db)2 

CONLC  and  CONRC  are  calibration  constants  (db)  obtained  from 
Calibration  Record  Words  624  (NB  LC),  625  (NB  RC), 
627  (WB  LC),  and  628  (WB  RC) 

POWERT  (peak  transmit  power  in  dbw)  for  NB 
=  PWRCN  +  PWRSN  log  XPKPWR 

POWERT  for  WB 

L  =  PWRSN  +  PWRSW  log  XPKPWR 

where 


PWRCN  is  Calibration  Record  Word  620 
PWRSN  is  Calibration  Record  Word  621 
PWRCW  is  Calibration  Record  Word  622 
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TABLE  I 


AVERAGING  PARAMETERS  FOR  ALC10 
Pulse  No.  of  Pulses  Averaged 


PRF 

Spacing 

u 

0.01 

u 

0.02 

u 

0.04 

u 

0.05 

u 

0.08 

o.iff 

(PPS) 

(s/pulse) 

200 

0.01## 

1 

2 

4 

5 

8 

10 

160 

uu 

0.0125 

(l)f 

(1) 

(1) 

4 

(1) 

8 

100 

0.01 

1 

2 

4 

5 

8 

10 

80 

0.0125 

(1) 

(i) 

(1) 

4 

(i) 

8 

50 

0.02 

(1) 

i 

2 

2-3n 

4 

5 

40 

0.025 

(1) 

(i) 

(1) 

2 

(1) 

4 

25 

0.04 

(1) 

(i) 

1 

(1) 

2 

2-3tt 

20 

0.05 

(1) 

(i) 

(1) 

1 

(1) 

2 

Averaging  interval(s)  (AVINT). 

## 

When  ALCOR  transmits  either  200  or  160  pps,  the  ADT  tape  (prior  to  April  71)  has 
always  contained  every  other  pulse  for  a  recorded  PRF  of  100  or  80  pps. 

+ 

If  AVINT  is  not  obtainable  at  the  PRF  without  using  fractional  pulses,  only  one  pulse 
is  averaged  (AVINT  =  1/PRF). 

++ 

For  this  AVINT  and  PRF,  2  pulses  and  3  pulses  are  averaged  alternately. 
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PWRSW  is  Calibration  Record  Word  623 
XPKPWR  is  ADT  Record  Byte  >44 

R,  Az,  and  El  are  corrected  as  follower 

R  =  IRANGE  +  TRBIAS  +  TTCOR  +  RRCOR  -  RCORF 

Az  =  IAZ  +  AZBIAS 

El  =  IEL+  ELBIAS  -  ECORF 

where 

IRANGE  is  uncorrected  R 

TRBIAS  is  range  bias 

TTCOR  (transit  time  correction)  =  RR/c 

RRCOR  is  range  doppler  coupling  correction 

RCORF  is  tropospheric  refraction  correction 

IAZ  is  Az  encoder  angle 

AZBIAS  is  Az  bias  (Calibration  Record  Word  602) 

IEL  is  El  encoder  angle 

ELBIAS  is  El  bias  (Calibration  Record  Word  603) 

ECORF  is  tropospheric  refraction  correction 

Alt  is  computed  as  follows: 

Alt  =  (R2  +  R  2  +  2  RR  sin  El)'*  -  R 
e  c  e 

where 

Rg  =  radius  of  earth  (6378. 145  km) 

Before  processing,  the  main  program  checks  that  ITBAND  (tape)  =  IBAND 
(input).  This  determines  that  if  WB  data  is  requested,  WB  data  exists  on  the  tape 
requested. 


III. 


OPERATION 


A.  Input 

Number  of  c.a?  s 

Launch  time  (total  GMT  ms) 

Range  offset  (in) 

Waveform  and  polarization 
Averaging  interval  (s) 

Start  and  stop  range  gates  for  peak  RCS  search 

First  and  last  pulse  nos.  of  processing  intervals  and  initial  PRF 

No.  of  processing  intervals 
A  sample  input  is  shown  in  Appendix  A. 

CARD  1  (12) 

(Col.) 

1-2  ICASEN  No.  of  cases 

CARD  2  (110,  F10.1,  515,  F5.2,  IX,  A4) 


1-10 

ILNCH 

Launch  time  in  total  GMT  ms 

11-20 

XlNOFF# 

Range  offset  (m)  (0) 

21-25 

NCELL1# 

Initial  gate  for  peak  RCS  plot  (50) 

26-30 

jf 

NCELL2 

Final  gate  for  peak  RCS  plot  (54) 

31-35 

NBAND 

0  =  NB;  1  =  WB 

36-40 

IPO  LA  R 

0  =  LC ;  1  =  RC 

41-45 

NVAIJS 

No.  of  processing  intervals 

46-50 

A  VINT  # 

Averaging  interval  (s)  (0.05) 

52-55 

TITL 

Title  for  listing  and  plots 

(U)#  If  left  blank,  program  sets  to  indicated  value. 
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CARD  3  (6110) 
(Col. ) 


1-10  NSTART  First  pulse  no.  of  initial  processing  interval 

11-20  NSTOP  Last  pulse  no.  of  initial  processing  interval 

21-30  INPRF  PRF  at  first  pulse  of  initial  processing  interval 

31-40  NSTART  First  pulse  no.  of  second  processing  interval 

41-50  NSTOP  Last  pulse  no.  of  second  processing  interval 

51-60  INPRF  PRF  at  first  pulse  of  second  processing  interval 

Repeat  Card  3  as  necessary. 

Repeat  Cards  2  and  3  for  each  case. 

B.  Output 

LISTING 


Selected  input  parameters 
No.  of  points  to  be  plotted 
Processing  interval  (TAL) 

PLOTS 

RCS  vs  TAL  and  Alt 


Sample  outputs  are  given  in  Appendix  B. 


IV. 


PROGRAM  LIMITATIONS 


NVALS 

XIVALS 

IPRS 

XINOFF 

No.  of  averaging 
intervals 


^  50  processing  intervals 

Either  0.01,  0.02,  0.04,  0.05,  0.08,  or  0. 1  s 

either  200,  160,  100,  80,  50,  40,  25,  or  20  pps 
(any  other  PRF  will  give  single  pulse  averaging) 

Must  be  on  tape 

s  6000 
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V. 


PROGRAMMING 


A.  L10ALC  (see  Appendices  C  and  D. ) 

L10ALC  is  the  control  section  of  ALC10.  L10ALC  reads  the  input 

C  ardS’  Calls  READGS  and  UNPACK,  and  averages  the  data  returned.  L10ALC  also 
calls  the  subroutines  that  plot  and  print  the  data. 

B-  HEDADT  (see  Appendix  E. ) 

Subroutine  HEDADT  unpacks  the  ADT  header  record  which  contains 
ndwidth,  reel  no.  ,  WTR  no . ,  date  of  mission,  and  mission  designator.  The  call 
sMtement  is  HEDADT  [  ISIG  ,  INBUF(l),  IEQM(l)  ] 


C. 


INPUT 

INBUF(l)  First  word  in  the  ADT  header  record  ## 

OUTPUT 


IEQM(l) 

IZBAND 

(bandwidth:  1  =  WB: 

IEQM(2) 

ITREEL 

(reel  no. ) 

IEQM(3) 

IT  WTR 

(WTR  no.) 

IEQM(4) 

IMTH 

IEQM(5) 

IDAY 

(Date  of  test) 

IEQM(6) 

IYR 

IEQM(7-9) 

ITDESG 

(mission  designator) 

READGS  t 

(see  Appendix  F. ) 

TKe  first  call  to  subroutine  READGS  opens  the  file  and  reads  the  ADT 
header  record.  The  second  cail  to  READGS  reads  the  ADT  calibration  record  and 
stores  the  values  in  a  buffer  area.  L10ALC  extracts  the  individual  calibration  values 

“  requires.  Each  subsequent  call  to  READGS  reads  an  ADT  data  record  consisting  of 
eight  ALCOR  pulses. 


Not  used. 

## 

^  INBUF  (2)  to  INBUF  (1803)  contain  the  remaining  words  in  the  record 

Same  as  READJS^  except  rewind  procedure  is  included. 
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D  UNPACK2 

Subroutine  UNPACK  unpacks  the  raw  data  from  the  ADT,  and 
translates  it  into  a  format  usable  by  the  IBM  360/67  computer. 

E.  REFC  (see  Appendix  G.) 

The  tropospheric  refraction  correction  subroutine,  REFC,  is  based 
on  tropospheric  refraction  tables  in  PPP-36J.  A  modified  version  of  this  subroutine 
is  now  in  use. 


The  call  statement  is  REFC  (E,  R,  DEE,  DRR) 

E  Uncorrected  El  (must  be  between  0°  and  90° 

R  Uncorrected  R  (ft) 

DEE  El  tropospheric  correction 
ERR  R  tropospheric  correction  (ft) 

The  corrected  values  to  be  computed  after  exiting  from  the  REFC 

routine  are: 


El  =  E  -  DEE 
R  (ft)  =  R  -  DRR 

F.  CELPLT 

Subroutine  CELPLT  plots  the  RCS  vs  TAL  plots. 

G.  Plotting  System  Subroutines 

They  are  REREAD,  STOIDV,  and  PLTND. 
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APPENDIX  A 
ALC10  INPUT 


1 


CARD  1 


0  0  0  0  0  0  0  0  0  C  0  0  0  0  0  0  0  0  0  0  0  0  0  0  B  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  D  0  0  0  0  ",  0  0  0  0  0  0  0  D  0  0  0  0  0  0  0  0  0  0  0  0  0 1)  0  0  0  0  0  0  0  0  0  0 
1  !  1  '  !  1  '  1  < 11 "  i!ijui5iiiniiimi!!aaK:i!imi»ii  usuim  min  me  «u>o  iiisiii.ihi  si  S'USi  si  is  si  srsisi  mi  mu  ms  ii  si  ms  Jin  hi) ms n  inmio 

1  I  I  I  1  I  I  HI  I  I  1  I  I  I  I  1  I  I  I  I  I  I  I  I  I  1  I  III  II  I  I  I  I  I  I  I  I  I  1  I  I  I  1  I  I  11  I  I  111  I  I  I  I  i  I  I  1  |  |  |  |  |  |  1  ]  1  |  |  |  |  | 

2  2  2  2  2  2  2  2  2  i  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2 2 2  2  2  2  2  2 2  2  2  2  2  2  2 2  2  2  2 2  2  2  2  2  2  2  2  2  2  2 U 2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2 

33333333333333333333333333333333333333333333333333333333333333333333333333333333 
<4<44444MM444M44<44444444H44444444444J444M4444<444444444M4MM4MMM44M4 
5  5  5  5  5  5  5  5  5  5  5  5  5  $  5  5  5  5  5  5  5  5  S  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  S  S  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 

fi  6  6  6  6  6  6  6  6  6  i  6  6  6  E  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  5  6  8  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  8  6  6  6  6  6  6  6  6  6  6  6  6  G  6  6 

77777777777/77777777777777J777777777777777777777777777T7777777777777777777777777 

8  8  R 8  8  8  8  8  8  C  8  8  B  8  8  8  8  0  6  8  8  8  ! !  8  4 8  1 4 ii  0 8  S  i  I 8  1 8 ! 8 0 8  8  3  8  9  9  9  0 6 ! 8  B  B  B  B  6  9  4  8  8  B  4 9  B  8  B  8  8  8  0  B  8  B I 0  8  B  8 B 
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APPENDIX  C 

L10ALC  PROGRAM  LISTING 

DOUBLE  PRECISION  XLNCH, 01000 
C 

DIMENSION  NSTART(50),NST0P(50> 

DIMENSION  XLCSUMt 170) ,XRCSUM( 170 ) , I NPRF  I  50  I  , 
l  XATBL 1 128) .XLCDB 1 170) ,XRCDB( 1701 , I LC AMP < 170  ) , I RC AMP ( 170) 

2,  ILCPHA (  170  )  ,  IRCPHA ( 170  I 

DIMENSION  IPRS(8)*101(8,2),I02(8,2)»I04(B*2)»I05(B»2)*IOB(B,2), 

I  110(8,2>,IUSE(B,2),XIVALS(6) 

DIMENSION  XNHUF I  1603) ,PIFA( 16) »OIFA( 16)  ,  XKRCSI5) 

DIMENSION  I ECM ( 9 ) , 1 TDESG ( 3  I 
DOUBLE  PRECISION  TCOM 

COMMON/ I  COM/ INPUF (1B03I,IAZ,IEL,INDEX,IPPRCS,I0RS, I  RANGE, I PKPWR , IR 
1D0T , I  ALT , INDAZ, JNDAZ, INDEL,IRB54,IRBR5,I0PRCS, I  240B 1, l 24CB2  * I  2  40B  3 
1, l 24  IB  1  *  I241B2»1241B3,XPPAGC, I  BE  T  A , ME  W A , I  BAND, NSW , R8 I  AS ( B 1 , I SVPR I , 

II  HRS, I  MIN, I  SEC, IMS EC, ISTAT (21 I.TRBIAS, ISTAT1*ISTAT2,ISTAT3,ISTAT4, 
1IALSW, ISTSW.NBWE, ISIGNO, 127812, JCON , NB EG , NEND , I T ST , NLMPR I ,XCPAGC, 

I  IT  BAND, I TAPNO, IPRF, I  POLAR, ISSERR,PIFA,CIFA,PFSA,OFSA,PSSA,CSSA, 
1PSSL.0SSL , ICCDF, I273B5, I273B6, 127387, I 273B8 , I MCVP , I MCVC , I OTFST , 

COMMON /PLOW /AL  SAV( 100 ) ,OBT  ,QBB  ,IRUN,K0UNT,NALT*NEWREC*PL0T10(B 
1) .CSRV(dOOO)  ,  KOM  <6000  ),TPINCH,TSAV(  IOC),  DP  INCH, ICASE 

EQUIVALENCE  (ILCAMP(1),IDAT<1))*{ILCPHA(1),IDAT(171)),(IRCAMP(1), 

II  DAT ( 341  )  )  , (  IRCPHAI  1)  ,  l OAT ( 51 1)1 
EQUIVALENCE (XNPUFI 1 ), INBUFI 1 II 

EQUIVALENCE  (  IFCM (  1  ), I  ZB AND) , ( I EQM (2),ITREEL),(IEQM(3),ITWTR), 

2 (  I  EQM ( 4  )  , I M T H  )  ,  (  I  EQM ( 5 ) , I D AY  )  , (  I  ECM I  6  )  ,  I YR ) , 

3 (  I  ECM ( 7 ), I TDESG I  1 1) 

DATA  ZLC/'LC  */,ZRC/'RC  '/.ZWB/'WB  */,ZNB/»NB  •/, BLANK/'  '/ 
DATA  I  PR S/ 200  i  16C, ICO, 80, 50, AO, 25, 20/ 

DATA  X  IVALS/.OI , . 02, . 04 , . 05 , . 08  ,  .  10/ 

DATA  101/  1,  1,  1,  1,  1,  l,  l*  1*  1*  l*  l»  l»  l»  l/ 

DATA  102/  2,  1,  2,  1,  1,  1,  1,  I,  2»  l*  2,  1,  1,  1,  1,  1  / 

DATA  104/  4,  1,  4,  1,  2,  1,  1,  l,  4,  1,  4,  1,  2,  1,  1,  1/ 

DATA  105/  5,  4,  5,  4,  2,  2,  l,  l,  5,  4,  5,  4,  3,  2,  1,  1/ 

DATA  108/  8,  1,  8,  l,  4,  1,  2,  l,  8,  l,  B,  1,  4,  1,  2,  1/ 

DATA  I10/1C,  8,10,  8,  5,  4,  2,  2,10,  B,10,  8,  5,  4,  3,  2/ 

DATA  ER  /6378.145/ 

DATA  D1000/1C00.  DO/ 

I  POL  AR  =  0  LEFT  CIRCULAR  DATA  REQUESTED 
I PCJL  AR  =  1  SIGHT  CIRCULAR  DATA  REQUESTED 
NB  AND  =>  0  NARROW  BAND  DATA  REQUESTED 
NBANO  =  1  WICE  BAND  DATA  REQUESTED 

NEW4  =  0  M I S  S  ION  F.OWN  BEFORE  15  CCT  70  (OLD  ATTN.) 

NEWA  =  1  MISSION  FLOWN  AFTER  15  CCT  70  (NEW  ATTN.) 
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CALL  STOIOVI * AIC1C60’ ,6,C ) 

CALL  RfcREAD  (9<H  530  ) 

K0UNT=0 
[ C ASE=0 

READ! 5  *  1 3 ) ICASFN 

13  FORMAT ( 12) 

[  F (  ICASEN.LE.O) ICASEN*1 
5  CONTINUE 
IFRST 1=0 
[FRST2=0 
[  FRST3  =  0 
[ F  RST4  =  0 
[  A  V  =  0 
JJ  =  0 

[CASE* [CASE* 1 

[FI  [CASE.GT. ICASENJGO  10  127 

REACH  5  *  1 )  [LNCH,XINGFF,NCELL1 ,  NCELL2 , NBAND ,  I HL AR , N VALS  ,  AV  I  NT  , 
ITIlLt  (NSTARTm  tNSTOP*  *),  INPRFI  1) ,  1*1,  AVALS  ) 

1  FORMAT  (I  [t0,Fl0.i,5[5,r>.2.1X,A4/(6ll0)) 

I F ( AV INT.Lr.O.  ) A  V  I N  T= . 05 
IE0F--0 
[  ERR=0 

CALL  REACGSI INRUF, [EOr, [ERR) 

I F ( ICASE.GT.l )00  TO  14 
IS [0*1 

CALL  HEOADT  (  I S I G , I NBUF ( 1 ) , I E CM (  1 )  ) 

[ T  BAND* I ZBANC 
NE'WA=0 

IF ( [VR.GT. 70 )GT  TO  282 
I F ( [YR.LT,70)GC  TO  283 
IF  (  [MTH.GT.IOGO  TO  282 
lr ( IMTh.LT . 10 )G0  TO  283 
[FI 10AY.LT. 15)00  TO  283 

282  NEHA  = 1 

283  CONTINUE 

14  I E  RR  =  0 

CALL  REACGSI INRUF, [EOF, IEKR) 

[FI ICASE.GT . 1  )G0  TO  12 

STORE  THE  DFS1RE0  CALIBRATION  VALUES 

N  =  0 

DO  20  K  =  256 »  383 
N=N*1 

20  XATBON)  =  XNBUF(K) 

N=  0 

DO  22  K=512,527 
N=N*1 

22  PI  FAIN) -XITBUF(K) 

N  =  0 

DO  53  K =528,543 
N=A  ♦  l 

53  OlFA(N)=XNeUF(K) 

PFSA=XNBUF(592) 
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PSSA  =  XNBUF<  5931 
0FSA=XNBUF(595) 

0SSA=XNBUF(595) 

AQIAS*XNBUF(60?) 

EBIASxXNBUF(607) 

OEGCON* »  180. *.05793691/3151. 59 
AZBI  ASxDEGCON*ABl AS 
ELPIAS= JEGCQN*FBI AS 


N  =  0 

DO  55  >'=605,611 
NXN  +  1 

55  R8  I  AS ( N I  = INBUFIKi 

PWRCN=XNBUF(620) 

PWRSN=XNBUF(621  ) 

PWRCW=XN8UF(62?) 

PWRSW  =  XNBU<:I623) 

N  =  0 

DO  57  K=625,62P 
N>N+  1 

57  XKRCS(N)xXNBUF(K) 

PSSL=XNBUF( 629» 

OSSLxXNBUF(630) 

12  LGOxO 

DO  8  K=  1 « 6 

IF ( ( XI VALSIK l-AVINT )  .L E . . 005  I L GO*K 

8  CONTINUE 

IF f LGO.NE.OIGO  TC  9 
WRITE(6,793»AVINT 

793  FORMAT ( •OILLEGAL  AVERAGING  I N T FR VAL ' , F  1 2. 3 ) 
GO  TO  127 

9  GO  TO  ( 21,23,25, 27,29,51) ,LGO 

21  DO  31  N= l »  2 
DO  31  M= l i 8 
31  lUSfcIH-N)xIOHM,N) 

GO  TO  52 

23  DO  33  N=  1 » 2 
DO  33  M= l »  8 
33  IUSE(M,N)X102(M,M 
GO  TO  52 

25  DO  35  N=  1  *  2 
DO  35  M= 1 »  8 
25  IUSE(M»N)X105(M,NI 
GO  VO  52 

27  DO  37  N*l,2 
DO  37  M*l»8 
37  IUSE(M,NI*I05(M,N) 


GO 

TO 

42 

29 

DC 

39 

N=  1 1 

,2 

00 

39 

H*  1 , 

>6 

39 

IUSE ( H 

liN)*I08(M 

GO 

TO 

42 

51 

DO 

41 

N  =  1 1 

,2 

DO 

41 

M=  1 , 

, a 

41  I  USE ( M  t N ) *  I 1C  I  M ,  N  ) 

GO  TO  42 

C 

42  CONTINUE 
C 

JCCN=-1 
I NDEX  =  0 
1 T  ST*  1 
I  PUL  S  =  0 
C 

DO  120  IJ*1,NVALS 
NBEG=NSTART( I J ) 

C 

DO  7  K  =  1 » 170 
XLCSUM ( K  )  *0. 

7  xpcsumixi=o. 

I PRF*  I NPRF ( IJ) 

I  N  T  A  V  =  1 

I F (  ABS(XIN0FF».GT.I)GC  TO  61 
OU  8S  N*l,8 

IF!  IPRS(N).NE.TPRF) GO  TO  85 

intavkusein.ii 

85  CONTINUE 
81  JX=1 
GO  TO  3 

2  CONTINUE 
JX  =  JX+  1 

I NTAV=  1 

I F (  A«S(XIN0FF).GT.1)G0  TO  3 
OU  83  N*l,8 

IK IPRSINI.NE.IPRFIGO  TO  83 
IN  TA V=  IUSE ( Nf  JX ) 

I F ( JX.EQ.2)JX=0 
83  CONTINUE 
C 

3  JCCN* JCON*  I 

»F(  JCON.EQ.9.QR.JCON.EC.OGO  TO  97 
INCEX* ( JCON-  11*900 
GC  TO  99 
97  JCCN=  1 
INC* X=0 
96  I  EOF  *0 
I ERR*0 

CALL  REACGSI INPUF, lEOFf IERRI 
I F ( I  EOF .NE .0 )G0  TO  127 
IF(  IERP.  .EC.1IGP  TO  103 
99  CALL  UNPACK 

XOFFST=( FLOAT  I IOF F S T ) /204  6 . 1*14.969625 
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JOFFST=XOFFST 

3 F  (  ICOOE.EC.5)YOPAGC*XPPAGC 
IF  IIC0CE.EQ.7)X0PAGC*XPPAGC 
101  IF  (  IFRST  2. EQ. 1 )G0  TO  92 
Z  BAN  =  ZN  B 

IF( ITBAND.EQ. I )ZBAN*ZWB 
ZPOL=ZLC 

l  F (  IP0LAR.EQ.1)ZP0L*ZRC 

IF(NCELL1.LE.0)NCELL1*50 

IF(NCELL2.LE.0)NCELL2*54 

PWRUS 1=PWRCN 

IF ( NBAND.EQ. I )  PWRUS 1*PWRCW 
PWRUS2=PWRSN 

IFI NBAND.EC. 1 )  PWRUS2*PWRSW 
CCNLC*XKRC  S (  II 
CCNROXKRCSI2) 

IFI NBAND.NE. 1 1^0  TC  17 
CCNLOXKRCSI  A) 

CCNRC  =  XKP.CS(  5) 

17  CONTINUE 

IFI ICCOE.EC.51CONRC=CONLC 
IFI IC0DE.EQ.7JC0NRC*C0NLC 

WRITE  I  6, 200  )  Z  PPL.  Z  BAN.  I TREEL , T I TL »  1 1 EQR II) , 1*4 , 6) 

200  FORMAT  I • 1L60-ALC0R  POLAR  *  • , A2 , AX, *  BAND  *  * ,A2,4X, ’REEL  NG.  *  * 
1,15,’  TITLE  =  • , A A ,  •  DATE  *  • , 1 2 , • / ' , 12 » ' / * » 1 2 ) 

IFI  ICASE.GT.lIf.O  TO  207 

207  WRITE(fc,208)NCFLLl,C0NLC, NBANO 

208  FORMAT  I • OCELL  REG  IN  I  SEARCH  )*  *,I3,5X 
l'KRCSILCI  *  •  , F  1C. 3,2X, • IBAND  *  *,Il,5X) 

WR I TE I  6  »  210  I  NCELL2  tCONRC ,  IPPIAR.NVALS 

210  FORMAT  I ’OCELL  END! SEARCH  I  )  *  *,I3,5X 

l’KRCS(RC)  *  *,F10.3,2X>  UPJIAR*  • ,  1 1 , 5X  ,  •  NV  ALS  *  *,I2) 

WR ITEI 6,2121  I NSTART I  I  I .NSTOPI II , INPRF(I) , 1*1,  NVALSI 
212  FORMAT ( ' 0  START  STOP  PRF 12X ,* START  STOP  PRF • ,  1 2X, ' ST ART 
1ST0P  PRF ', 12X, ' START  STOP  PRF*  / 1  A I 2X , I  5 , 2X , 1 5 , 2X , I  3, 10X ) ) ) 
XLNCH*DFLOAT (  ILNCHI/01000 
WRITEI6, 214 )XLNCH,AVINT ,XINOFF 

2 1 A  FORMAT  I ' 0  LIFT  OFF  i 1  ME  *  ’,F1C.3»’  AVERAGING  INTERVAL  (SE 

ICS)  *  ’,F8.2»*  INPUT  RANGE  OFFSET  !W)  *  SFIO.IJI 

RRUSE*- . 0094  3 

IFI ITBAND.EG*l )RRUSE=-. 000115 
IFIN8AND.NE. ITPANDtGO  TU  695 
I FRST  2=  l 
92  CONTINUE 

620  IF  I NUMPRI.LT. NSTART I  I J ) )G0  TO  3 

IFI  ABSI XINOFF -XCFF ST ) . GT . 1 . ) GO  TO  117 
I TST= 1 

627  IPULS* I PULS+1 

IF  I  IPULS.LT. INTAVIGO  TO  10 
I  FULS  =  0 
I TST*2 

10  CONTINUE 

IFI  IFRSTA.EQ. I )GC  TO  610 
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IPRCLD=  IPRF 
IFRST4= l 

GO  TO  61  l 

6 1 C  !FI  IPRF.EC.IPRPLOGO  TC  611 

WR!TE(6,62A1NUMPRI,  IPROLD.IPRF 

62 A  FORMAT ( /25X’CU°RFNT  PRI  «  *,I8  .*  OLD  PRF  *  tI5*  CURRENT 
1  ’,15/) 

611  CONTINUE 

IK  I SSERR.NE.  1 1GC  TO  617 

612  FORMATI  ^SLOh^HITCH  BITS  ARE  BOTH  *  Q*,5X,*  CURRENT  PRI  * 

617  CONTINUC 

IT0T=(3600*IHRC  +  60«'IMIN*1SEC)*1000*IMSEC 

I  r AL  = I TOT-ILNCH 
TAL=DFLOAT( ITAL )/DlOCO 
TOTL=DFLCAT(  ITPTI/DIOOO 
RDCT*(  IRDOT/I 8192.0) )*14.  989625 

RANGE=  (FLOAT  (  l<>  AN  GE)/ 2040000.  )* 14 . 989625* TR B I  A S* .  14989625 
T  T  CCR  r  (RANGE /2997T6.  )  * ( RDO  T /1C0C. ) 

RANGE=RANGF*T  TCOR 
RRCCR  =  RRUSE*RDr'T 
RANGE=RANGF+RRPOR/ 1C00. 

E  L  =  ( IEL*2*3141.59266)/(2. 0**17) 

XEL=EL*.C57295P 
XFL=XEL*ELBI AS 

CALL  REFC(XEL,OANGE,ECCRF,RC0RF) 

RNGF=RANGE-RCO°F 

elvf=xel-eccrf 
R  AOEL  =  EL  VF  * • 0 1  7453 

Al T  =  SCR  T (  RNCF**2*ER*ER*2. *RNGF  *cR  *  S I N ( R ADEL  )  )-ER 

range=rngf 

XT  HR  =  4C . * ALUG  lr ( RANGE  ) 


X  PKPWR  = IPKPWR 

P0WERT=PViRUSl*t’WRUS2*AL0G10(X  PKPWR) 


C 

IF ( IFRST l.EU.l >GC  TC  11 
talueg=tal 

NPRBEG=NUMPR I 
AL  TBEG=  Al T 
IFRST  1=1 
ll  CONTINUE 
C 

I  A V= I  A  V  *  1 

I F ( I  POL AR.EC. I  1 GC  TO  73 
DO  79  K=NCELL1,NCELL2 

N*  I  LCAMP (  K  )  *  1  _  . 

XLCDB ( K ) *X ATBL ( N 1 ♦ XTRR* XPP AGC+CONLC-POKERT 

XLCSUM(K)=XLCSl'H(K)*10.**(XLCD6lK)/10. ) 

79  CCNTINUt 
GO  TC  74 

73  DO  76  K-NCELL1.NCELL2 

M=  IRCAMP (  K  )  ♦  1  „  _ 

XRCDB(K)*XATBL  ( H ) *X TR R* XCP AGC+CONRC-POWER T 

XRCSUH(K)*XRCSUM(K)*10.**(XRCDB(K)/10. ) 

76  CONTINUE 
C 

74  GO  TO  (  199,  19) , I TST 


PRF  « 


’,110) 
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199  l  F  (  NUMPR  I-NE-NSTOP(  IJMGO  TC  117 

19  CCNTINUE 

XDBPK*- 1000 • 

IF(  IPOLAR.EQ. 1 )G0  TO  77 
00  72  K=NCELL 1 , NCELL2 

XLCCB!K)=10.*AL0G10(XLCSUM(K)  /FLOAT ( I AV  )I 
IF (XLCDB(K).LT . XDBPK ) GO  TO  72 
XDBPK*XLCDB ( K  ) 

I RGPK*K 
72  CCNTINUE 
GO  TO  93 

77  DO  78  K-NCELL1,NCELL2 
XRCDB!K)=10.*AL0G10(XRCSUM{K)  /FLOAT! IAV  )) 
IF(XRCCB(K).LT .XDBPK) GO  TO  78 

XDBPK  =  XRCDB ( K  ) 

I PGPK*K 

7 8  CCNTINUE 
93  JJ=JJ+l 

I RUN* J J 

25b  IF ( IRUN.GT. L )GC  TO  260 
NALT  *  1 

T  S  AV  (  N  AL  T  )  =  T  AL 
ALSAViNALT)*ALT 
GO  TO  360 

260  IALT1=  IFIX(ALPqtV  /  10« ) 

IALT2=IFIXIALT/1C.) 

IF!  IALTI.EC.IALT2JGC  TO  360 
IF! ALT .GT , ALPRFV ) GO  TO  280 
A I.  TEST* I  ALT  1*10 
GO  TO  300 

280  AL7EST*  I  A‘..T2*10 
300  N AL  T  =  NAL  T ♦ 1 

ALSAV(NALT) =  AL  TE  ST 

FAC=! ALTEST-ALPREV  I/! ALT-ALPREV) 

TSAV!NALT ) *T PRF V* F AC* 1 T AL-TPRE V ) 

360  ALPREV=ALT 
TPREV  =  T  AL 
TCOM! JJ )=TALBEO 
CSRV! JJ)=XCBPK 
71  00  82  K* l » l 70 
XLCSUM ( K ) *0 . 

82  XRCSUN ! K ) *0. 

I A  V  =  0 
IFRSTA=0 
Ir-RSTl  =  0 

IF ( JJ.GE.60C0 )G0  TO  121 

118  IFINUMPRI.LT. NSTOP! IJ))GC  TC  2 
117  IFINUMPRI.LT. NST0P!1J))GC  TO  3 

119  I FRST  3  =  0 
I  PULS*0 

I A  V  =  Q 
IFRST4=0 
I FRST  1  =  0 

120  CCNTINUE 

121  CONTINUE 
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NALT *N  AL  T  l 
TSAV!NALTJ=TAL 
ALSAV!NALT  ) =  AL  T 
NEWREG=0 

GG  TO  125 

103  WRITE!  6, 1071NUMPRI 

107  FORMAT ( 'OPARITY  ERROR  ON  READ  AFTER  PRl  *  * , I  I C ) 

GO  TO  9 <3 

680  WRITEI6, 109 ) NUMPR  I 

109  FORMAT!'  END  OF  FILE  REACHED  LAST  NUHPRI  VALUE  *  ',110) 

GO  TO  125 

695  WR  ITE! 6,  114 ) I8AN0, l TBANO 

114  FORMAT!'  INPUT  BAND*  ' I  1 0 »  *  BAND  ON  TAPE  *  *110) 

GO  TO  127 

125  PLOTID! I ) * T l TL 
PLOT  I  0 ! 2 ) =8L  AN* 

PLOT  10(3) =2  BAN 
PL0TID!4)*ZPCL 
00  126  K*5 , 8 

126  PL0TI0!K)=BLANK 

T  M I N=  T  COM  l  l  I 
TMAX  =  TCOM  (  ,1 J  ) 

WRITE <6,66 ) I  RUN, TMIN,  TMAX, I C ASE , NC ELL l , NCELL2 
66  FORMAT!  'OPLOT',15,'  POINTS  FROM  *,F15.3,'  TC  *,F15.3, 

11  OX, 'CASE  N3',T5,'  FROM  CELL  • ,  1 5  ,  '  TCCELL'cI5) 

640  CALL  CELPLT  TCDM VC SRV , J J , PLOT  1 0 , NCELL 1 , NCELL2 , ALS AV  , N AL T , T SA V, 
IXINOFF) 

I  E0F  =  2 

CALL  READGS! IN»UF, IEOF,  IERR  ) 

GO  TO  5 

127  CALL  PLTND 
RETURN 

END 
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APPENDIX  D-2 


APPENDIX  D-3 


STORE- 

calibration 

RECO  R.D 


/legal\ 

AVERAGING 

^INTERVAL- 


GENERATE 
PROPER 
PRF  TABLES 


< 


APPENDIX  D-5 
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PRF 

CHANGE 


PRINT 

^MESSAGE 

\  t 


APPENDIX  D-7 


INITIALISE 

PLOT 

TITLES 


* 

* 


XISIG 

XICAL 

XIEQU 

BASE 

HEOACT 


RETURN 


TEMPI 
TEMP2 
SAVE  A 


APPENDIX  E 

SUBROUTINE  HEDADT  PROGRAM  LISTING 

CALL  HEOADT  ( I S I C * INBUF ,  IE3U ) 

ISIG  *  1  UNPACK  THE  20  WORC  AOT  HEAOER 

START 

ENTRY  HEOAOT 
SPACE 
ECU  A 
EQU  5 
ECU  6 
ECU  12 
SPACE 

SAVE  flA.L2l.Tt* 

8 ALH  12,0 

USING  ♦.BASF 
ST  13.SAVEA+A 

LA  7 , SAVFA 

ST  7*8(0-131 

LR  13,7 

SPACE 

LH  XISIG, XIEGU, 0<  1) 

SPACE 

L  8,0(XICALJ 

ST  8, TEMPI 

ST  8.TEMP2 

SRL  8,31 

ST  8 , 0 ( X  I  ECU )  MBAND 

u  8, TEMPI 

GIL  8,1 

SRL  8,25 

ST  8 , A ( X l ECU  )  MREEL 

SPACE 

L  8 , A ( XICAL ) 

ST  8, TEMPI 

ST  8, TEMP2 

SRL  8,16 

ST  8 , 8 ( X  I  ECU )  MKTR 

L  8, TEMPI 

SLL  8,16 

SRL  8,2A 

ST  8,  1 2 ( X  I  ECU  I  MMNTH 

L  8.TEMP2 

SLL  8 , 2 A 

SRL  8,2A 

ST  6 ,  16  (  X  I  ECU  )  MOAY 

SPACE 

SR  8, 8 

IC  8 , 8 ( X I  CAL  I 

ST  8 , 2C ( X  I  ECU )  MYEAR 

MVC  2A(9,“IEGU)  ,9<XICAL>  MISSION  DES. 

SPACE 

L  13,  S AVEA  +  A 

RETURN  ( 1 A , I  2  ) ,  T 
CNUP  0  *  A 

DC  F  ’  0 ' 

DC  F  *  0  * 

DC  1 8 A ( *  I 

END 
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APPENDIX  F 

SUBROUTINE  READGS  PROGRAM  LF  1 1NG 

*  CALL  READGSt INBUF* IEOF, (ERR) 

START 

ENTRY  REACGS 
SPACE 

XZBUF  ECU  A 

XECF  EQU  5 

XERR  ECU  6 

BASE  EQU  12 

,PACE 

REAOGS  SAVE  (U, 12), T,* 

OALR  12,0 
USING  *,BASr 
ST  13 , SAVEA+A 

LA  7  *  SAVE A 
ST  7,8(0,13) 

LR  13,7 

SPACE 

LM  XZBUF, XERR.Ot 11 
L  7 , 0 ( XFOF  ) 

C  7, TWO 

BE  S9<3 

SPACE 

L  7, IFRST 1 

C  7 , TfcQ 

BE  I NI  T2 

C  7, ZERO 

BNE  WHICHF 

SPACE 

OPEN  UNDCP,  I  INPUT)  > 

INIT2  REAL  RCB3 , SF , INCCB, BUFFI, 7212 

^  HECK  R0B3 

HVC  BUFNUH  {  A  )  ,  ZERC 

HVC  IFRST1CA), CNE 

B  SKI 

SPACE 

WHICHF  L  3 ,  BUFNUH 

S  3, CNE 

BH  NEXTBUF  2 

B  NEXTBl'Fi 

SPACE 

NEXTBUF 1  HVC  BUFNUH ( A ) , ZERO 

CHECK  ROB  1 

SKI  READ  RDB2.SF,  INDCB,BUFF2, 7212 

L  q.AeUFFi 


NEXTBUH2 


LCCPC 


LOCP/ 


INDUP 


LUTLP 

BACRC 


ENCFILE 

S99 

RETURN 


INDCB 


ZERO 

CNE 

TWC 

IFRST  L 
BUFNUM 

ABUFF1 

ABUFF2 

SAVEA 

BUFFI 

t>UFF2 


B  l;opq 

SPACE 

M VC  DUFNUW (  4  ) , CNE 
CHECK  RCB 2 

READ  RCB1.SF, INDCB, BUFFI, 7212 

L  9 1  ABU^F  2 

SPACE 

LR  1C i XZPUF 

SR  11,11 

SR  3,3 

LA  8,180^ 

L  7 , 0  ( 3,9) 

ST  7,0(11,10) 

BCT  8, INDUP 

B  CUTLP 

LA  3,41 3  > 

LA  1 1 , 4  (  U) 

B  LCOPZ 

B  RETURN 

SPACE 

L  2, CNE 

ST  2.0JXFRR) 

BR  14 

SPACE 

L  2, ONE  STORE  END  OF  FILE  INDICATOR 

ST  2,0(XTCF) 

CLOSE  IINCCP, REREAD)  ,TYPE  =  T 
PVC  IFRSTK4) , TWO 

L  13,SAVEA*4 

RETURN  I  14,  12) , T 
SPACE 

DCB  CS0RG=PS,MACRF*(RC),DEVD*TA,DEN*2, BUFNC= 1,ECDAD=ENDFILE,C 
SYNAD=BADRD,BFTEK*S,DDNAME*FTllFC01 

SPACE 

CNOP 

DC  F  ' 

CC  F  *  l  * 

DC  F*2* 

CC  F ' 0 ’ 

DC  F ' 0  * 

SPACE 

CC  A ( BUFC  l ) 

DC  A( BUFF  2 ) 

EC  1 8  A ( ♦  ) 

SPACE 

CS  !.  8  0  3  F 

DS  1803F 

END 
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APPENDIX  G 


SUBROUTINE  REFC  PROGRAM  LISTING 


VERSION  6/16/70 


SUBROUTINE  KEFCSE.R.OEE.ORR) 

0  I  MENS  I  ON  OE(  16,8), CR<  16,8), EO(  16), RU<8) 

OUA  OE/O.C  ,0.0  ,0.0  .0.0  .0.0  .0.0  £.0 

MMBBB 

!  : ti-itssisiss&Kaa-.ss: 

E0.025A, 0.0121, 0.0  /  „  n  n  C.c,  0.0. 

OATA  OR/  C.c,  0.0.  0.0.  0.0.  0.0,  0.0,  O.U. 


1  0.0 

2  l  I  .A 


0.0.  ’■'.0. 

18.5,  IT. 6 


0.0.  0.0 

16.0.  16.1 


0.0, 

1 A  •  8  1 


C.O,  22.6,  21.5.  20. A,  19.9, 
IA.2.  13.2.  12.0.  1C. A,  a. 6, 


U:i: n.;;  .;:i; 

9. A, 


5  «:V,  i**T.  3*. a.  »»•*■  ii'V  ivy.  ia'.R',  iY.t,  26. A.  2a.o, 

6)40*0 1 !6?*C»l01»Ct  *  V  *  ’ ,  ,r  A.  06 • 3 »  T3.6«  56*®* 


l /. •  3  •  12* Tt 
) * • 4 i  31.4 


9. A,  8.1,132.0,  98.5.  T».A, 
26. A,  23.9,  2C.I.  16* A .  12.7. 


63.2, 

a.  1 1 


»  20.2 

a  96.  ) 


12.8, 

JT.7, 


9.5 

26.5 


16  A  l?  8.  9.5,  aii’.AOsIo’.irC.O.lOA.O,  86.3,  *3.6, 

8  A6.3.  \\W,  P‘.T.  26  5.  94.1.  20.3,  U.5.  ,2.8.  £». 

89*.a!  B*.a!aA6*.o',IT2. 0.103.0,  8T.A,  TA.O.  58.0,  A6.6, 

l0AT*JKC/’0°0U?r>!  ,3o!,60.  .200.  ,900.  ,1000.  ,2000./ 

IFIR.LE. 0.0)00  TC  300 
RG*R/l.852C*C0 
00  100  IE0-2.15 
{  «  l  f  -  I  E  U 

IFIE.GE.EOI  I  )  K-O  TC  1 20 
100  CONTINUE 

1*1 

120  00  200  JR0*2 , 8 

J.  10-JH0 

IFIRG.GE.RCU)  1GC  TO  220 
200  CONTINUE 

220  |r«alOG(RG/RCU))/ACOGIRO(J*1)/ROU)) 

I F ( E • L E . 0 .0 1  GO  TC  320 

0RR.(0R1*CK2*0°I i» J)) 

GO  TO  AOO 
300  0EE*0.0 

0RR=0 • 0 

320  OEE-OE 1  I , J)* IOC l I , J*1 )-UE 1  I » J>  > *2R 
ORK*OR(  I  ,  Jl  +  IORI  I  ..HU-OKI  !  .J)  )*2R 
GO  TO  AOO  , 

3A0  OELT*(E-EO<m/(EO(lMS  EO  I>> 

0EE*OELT*(DEI I* 1 » J)-OEt I , J )  )«0E< I , J ) 

ORR*OELT*(CR(  im.j)-orii.J)»torii.J) 

AOO  0RR*0RR*. 30A80-C3 
RETURN 
ENO 
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